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Abstract 

We consider the quasi-de Sitter geometry of the inflationary universe. We cal- 
culate the energy flux of the slowly rolling background scalar field through the 
quasi-de Sitter apparent horizon and set it equal to the change of the entropy 
(1/4 of the area) multiplied by the temperature, dE = TdS. Remarkably, 
this thermodynamic law reproduces the Friedmann equation for the rolling 
scalar field. The flux of the slowly rolling field through the horizon of the 
quasi-de Sitter geometry is similar to the accretion of a rolling scalar field 
onto a black hole, which we also analyze. Next we add inflaton fluctuations 
which generate scalar metric perturbations. Metric perturbations result in a 
variation of the area entropy. Again, the equation dE = TdS with fluctua- 
tions reproduces the linearized Einstein equations. In this picture as long as 
the Einstein equations hold, holography does not put limits on the quantum 
field theory during inflation. Due to the accumulating metric perturbations, 
the horizon area during inflation randomly wiggles with dispersion increasing 
with time. We discuss this in connection with the stochastic decsription of 
inflation. We also address the issue of the instability of inflaton fluctuations 
in the "hot tin can" picture of de Sitter horizon. 
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I. INTRODUCTION 



The inflationary paradigm established during the last 20 years assumes that the primor- 
dial equation of state is almost vacuum-like: p ~ — e. To realize this equation of state, most 
models deal with a scalar field <f>(t) (or other fields which in combination act as an effective 
scalar field) slowly rolling to the minimum of its potential V(<f>). During the slow roll regime 
the homogeneous scalar field produces geometry which can be well approximated by the 
quasi-de Sitter metric. 

The full pure de Sitter spacetime, which corresponds to a 4d hyperboloid of constant 
curvature, can be compactly represented by its Penrose diagram, given by the full square in 
Fig. |l|. It can be covered by different coordinates. Cosmologists most often use coordinates 
in which the metric is time-dependent and corresponds to an expanding flat universe 

ds 2 = -dt 2 + e 2Ht (dr 2 + r 2 dn 2 ) , (1) 

where dVt 2 = d9 2 +sin 2 9d(f> 2 . This coordinate system covers the upper half of the hyperboloid, 
which corresponds to the expansion branch. The Penrose diagram of de Sitter spacetime 
in flat FRW coordinates is shown on the left panel of Fig. [I]. Quasi-de Sitter geometry is 
described by the scale factor a(t) = aoe-f dtH ^\ where the Hubble parameter if is a slowly 
varying function of time, H <C H 2 . 

The time-dependent form of the metric ([I]) is very convenient for investigating the dy- 
namics of a scalar field with the equation 

□0 = V* (2) 

and for quantizing this field in the de Sitter spacetime . Among quantum scalar fields with 
mass m and conformal coupling £ in de Sitter geometry, the case of minimal coupling £ = 
and very small mass m <C H plays an especially important role. Indeed, the regularized 
vacuum expectation value is (Sep 2 ) = 2 . Formally, as was noted before the discovery of 
inflation, this is an odd case since its eigen-spectrum contains an infrared divergent term: 
{5(f) 2 ) —>■ oo as m —>■ 0. On the other hand, this is the most interesting case for application 
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FIG. 1. Penrose diagram of de Sitter spacetime in the flat FRW coordinates (left) and the 
static coordinates (right). Each point represent a sphere S 2 . Its radius at the horizon (dashed line 
on the left, edge of diamond on the right) is equal to jj. 
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to inflation, since the theory of inflaton (as well as tensor) fluctuations is reduced exactly 
to this case. Following the time evolution of individual fluctuations, it was found that the 
infrared divergence can be interpreted as the instability of quantum fluctuations of a very 
light scalar field, which are accumulated with time {S(f> 2 ) — §-t f|]|||4[]. Fluctuations of 
5(j) induce scalar metric perturbations 0|§||[7||| • This picture is a basis of the inflationary 
paradigm so successfully confirmed observationally. Notice that heavy or conformal fields 
are not produced by inflation. Further, backreaction of fluctuations 5(f) leads to the picture 
of stochastic evolution of quasi-de Sitter geometry PJlQl, and at large values of H even 
to self-reproduction (eternal) of the inflationary universe [IT]. Scalar field in the eternal 



inflationary universe is described naturally in terms of the probability distribution function 

p(0,t) eh. 

Recently, de Sitter spacetime and inflation have drawn significant attention in the theo- 
retical physics/superstring community. Some of the most interesting topics are holography 
and the thermodynamics associated with the de Sitter horizon. In this context, the static 
form of the metric of the de Sitter spacetime 

ds * = -(i- H * R 2 )dr 2 + {1 _^ 2R2) + R 2 dtf (3) 

is commonly used. The Penrose diagram of de Sitter spacetime in static coordinates is 
plotted on the right panel of Fig. |l|. The classical result of |13| is that observer at the origin 
detects a thermal radiation from the de Sitter horizon at R = jj with the temperature 
T = j-, and the horizon area A = jp is associated with the huge (geometrical) entropy 
S = ^7. Thermal vacuum in the causal patch ("hot tin can") corresponds to the Bunch- 
Davies vacuum of the metric (|]) 014,151 and gives a complementary picture of scalar field(s) 
fluctuations. It is not clear to us, however, how quantum fluctuations in the "hot tin can" 
picture correspond to the instability of quantum inflaton fluctuations 5(p and generation of 
metric perturbations. We will return to this point at the end of the paper. 

One of the issues in the holographic approach is the bookkeeping of entropy of de Sitter 
spacetime. The holography bound declares that the geometrical entropy of the horizon 
exceeds the entropy of quantum states (of fields and particles) within the volume surrounded 
by the horizon. It was recently claimed that counting the entropy of quantum fluctuations 
generated during inflation in the "hot tin can" and comparing it to the change of the apparent 
horizon entropy violates the holography bound unless an ultra-violet cutoff of order of ~ 10 16 
GeV in the momenta of fluctuations is imposed [JT5 . 



While it is expected that the approaches based on the time-dependent form of the de 
Sitter metric with unstable fluctuations and the static form of the de Sitter metric with 
thermal flux should give us complementary insights, their languages are apparently different. 
This is partly due to the difference between quasi-de Sitter and pure de Sitter geometries, and 
partly because different questions are addressed. However, we have to understand how these 
two different approaches to (quasi-) de Sitter geometry with a scalar field are compatible with 
each other with respect to such important issues as the generation of fluctuations, entropy 
and global geometry. 

In this paper we consider a particular question of how the apparent horizon area A, 
or the entropy S = vary due to the slow roll of the background scalar field and the 
generation of scalar metric perturbations during inflation. A novel element here is that we 
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combine the concepts of a dynamical, slowly rolling background field and the instability of 
its fluctuations, with the concept of geometrical, holographic entropy. 

In Section OL we will calculate a variation of the geometrical entropy due to the energy 
flux through the apparent horizon area. We find that, remarkably, the thermodynamical 
relation SE = TdS is equivalent to the Einstein equation for the rolling inflaton field. 
In a sense, our derivation of a correspondence between thermodynamics and the Einstein 
equations for inflation is a realization of such a correspondence found in an inspiring paper 



17] for local accelerating observers. However, we introduce a technique to treat the apparent 



horizon of R x S 2 topology which is different from the description [T7J of a local Rindler 
horizon for an accelerating observer. 

As we will see, a non- vanishing flux is generated by the kinetic term 2 of the slowly 
rolling inflaton field. It turns out that this problem is very similar to the problem of the 
interaction of a homogeneous rolling scalar field with a runaway potential V(<p) and a black 
hole. In Section |TTTj, we switch our attention from inflation to black holes. A rolling scalar 
field interacting with a black hole is a transparent illustration of the energy flow of a light 
scalar field through a horizon. 

In Section |TV|, we return to inflation. On top of the rolling background inflaton, we con- 
sider inflaton fluctuations 5<ft, which generate scalar metric perturbations <£>. We study the 
energy flux through the horizon including inhomogeneous 5<f) fluctuations and corresponding 
variations in the area of the horizon, or entropy dS, which are sensitive to the scalar metric 
perturbations In this case, the calculations are more involved than the calculations for 
the homogeneous time dependent background field in Section [H| This happens because 
there is no exact Killing vector generating the horizon. However, for metric perturbations 
which preserve spherical symmetry we still can define TdS and compare it with the energy 
flux through the horizon. 

In Section [V|, we argue that the metric perturbations generated from inflaton quantum 
fluctuations can indeed be treated as (locally) spherically symmetric. We apply the general 
formalism for spherically symmetric non-static geometry of Section |V| to fluctuations from 
inflation. Again, we find that the thermodynamical relation 5E = TdS leads to equations 
connecting $ and 5(j) which are in exact agreement with the linearized Einstein equations 
for the fluctuations from inflation. This allows us to give new insights into the entropy of 
cosmological fluctuations, as we will discuss in Section [V|. 



II. SLOWLY ROLLING INFLATON AND DE SITTER ENTROPY 

We begin with an inflaton scalar field in the spacetime ([I]). The homogeneous background 
field (p{t) is time-dependent and obeys the equation @, which for the metric ([[]) reads 
+ 3H(f) + = 0. If the Hubble parameter is large, the friction term here is significant and 
slowly rolls towards the minimum of V(4>), realizing the inflation. 

The apparent horizon of de Sitter spacetime has the topology of R 1 x S 2 , and is shown 
in Fig. [I] in both coordinate systems of interest. Imagine a spherical light front propagating 
from outside towards an observer at the origin r = or R = 0. In the time-dependent 
coordinates the space is exponentially expanding. There is a critical, trapped surface S 2 , 
which separates the region inside the apparent horizon which is in causal contact with the 
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observer, and the outside region, where expansion keeps the light from reaching the inside 
region. The instant physical radius of an S 2 sphere is equal to p = a(t)r. Since the position 
of the apparent horizon is r = , its radius ph is equal to jj. It is constant in the 
pure de Sitter space time, and a slightly increasing function of time in the quasi-de Sitter 
geometry where H(t) slowly decreases with time. 

So far we have not used the Einstein equations, which relate the energy momentum 
tensor = 0^0^ — {\4^^a — ^(0)) ^ °f t ne background inflaton field 0(i) to the scale 
factor a(t). Actually, we are going to derive them from the first law of thermodynamics of 
the de Sitter horizon. 

Let us calculate the flux of energy through the area A = of the apparent horizon of 
the quasi-de Sitter spacetime. In the spirit of the adiabatic approximation, we will treat 
the horizon as a static surface for this purpose, but allow it to slowly vary with time when 
calculating the change in its area. The energy flux through the horizon is given by the 
integral 

5E = 5jdE^C, (4) 

where effi^ is the 3-volume of the horizon, and £^ is the null generator of the horizon, as 
sketched in Fig. ^. Notice that at the horizon <iE M is parallel to £ M , as both are null. We 
will equate the energy flux through the horizon to the change of geometrical entropy dS 

TdS = 5E = 5 J T$ C , (5) 

where the temperature T is defined by the surface gravity at the horizon. For the de Sitter 
horizon, T — rr- 

The entropy dS and the expression (^) are geometrical invariants and can be calculated 
in arbitrary coordinates. In the case of de Sitter geometry we can work either in time- 
dependent (0) or static (|3|) coordinates. The metric ([I]) can easily be modified for the 
(adiabatic) regime of slowly rolling H(t), while it is not so straightforward for the static 
form of the metric. Therefore, we will work in time- dependent coordinates. Consider a 
slowly rolling inflaton 0(t) that is time-dependent but homogeneous in the coordinates ([[]). 
The background stress-energy tensor is diagonal. As the (approximate) Killing vector 
is £ M = (1, — Hr, 0, 0), the energy flux through the horizon T^^t^ = 2 is non-zero. To 
calculate the right hand side of equation (|^), we integrate over the infinitesimal volume 
47r dX, where A is the affine parameter along the null generators of the apparent horizon. 
In the left hand side of (|^), the variation of the entropy is dS = = — j^. By comparing 
both sides of the thermodynamic equation (^) and dividing them by dX = dt, we obtain 

H = -AtiG^) 2 . (6) 

This is nothing else but one of the Einstein equations. Using the equation of motion for 
and the equation (0), and assuming the field is dominant, one can reconstruct the 
Friedmann equation for H. Note that for a pure de Sitter geometry = and the energy 
flux through horizon is absent. 

Formula (§) can also be derived from ([|) in the coordinates (|3]). The coordinates (t,r) 
of (|ip and (r, R) of ([|) are related by the transformations 
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In static coordinates, we have <f> = <fi (r + 7^7 — H 2 R 2 )j. This 0(r, -R) dependence cor- 
responds to the field profile delayed as R approaches the horizon. We find a similar delayed 
field configuration in black hole spacetimes in the next Section. The energy flux through 
the horizon = iT 0,ij ~ (f 2 is non-vanishing. In these coordinates, however, a central 
observer does not see the flux going through the horizon but rather accumulating the energy 
just outside the horizon. To quantify this, one can define the local mass function M as it is 
done in Appendix B. 

To conclude this section, we consider a sub-dominant, test scalar field x during inflation. 
We can calculate the flux of the field x through the de Sitter horizon. In this case, it is 
not important whether we are dealing with pure de Sitter or quasi-de Sitter geometry. In- 
dependent of the potential V(x), the x-field contribution to the energy flux SE is given by 
precisely the same expression as the inflaton field contribution. Either in terms of the ther- 
modynamic equation TdS = SE or in terms of the Einstein equations (|j), the contribution 
of x is additive 

H = -4ttG(0 2 + X 2 ) ■ (8) 

This formula works even if the inflaton field velocity is sub-dominant, <fi <C x-> an d indepen- 
dent of the shape of the x-field potential. In particular, it works for a massive scalar field 
X, which oscillates about the minimum with an exponentially decreasing amplitude. 



III. ROLLING SCALAR FIELD AND BLACK HOLE 

In this section we interrupt the discussion of the slowly rolling inflaton in the early 
universe and focus on the seemingly different topic of a scalar field accreting into a black 
hole. The reader of the de Sitter story can jump to the next Section. 

Consider a homogeneous scalar field (f>{t) with a runaway type potential V((f>), say V ~ 
l/(p 2 or V((p) ~ e -<t>/ M p . The concrete form of the potential is not important. We can 
even broaden the class of models to include fields with the potential V(4>) = ^m 2 (fi 2 with 
very small masses m, as long as we deal only with the phase of slow roll towards the 
minimum. The well-known case of a massive scalar field oscillating around its minimum will 
have qualitatively different accretion onto a black hole similar to the accretion of massive 
(quasi)-particles. 

We will consider a test homogeneous (cosmological) scalar field described by the equation 
(§). In the vicinity of a black hole, the spacetime is described by the Schwarzschild metric 

, 9 / 2M\ , , dr 2 9 ,) , . 

ds 2 = - f 1 - — dt 2 + + r 2 dQ 2 , (9) 

V r ) 

and we will treat the cosmological evolution of the scalar field as the boundary conditions 
imposed sufficiently far away from the black hole. The scalar field retains its spherical 
symmetry = 0(t,r), so near a black hole it obeys the equation @) 
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Q r 2 



dt 2 



1 



2M\ 2M 



1 



2M 



V 6 = 0, 



(10) 



where the tortoise coordinate r* is defined as r* = r + 2Mln(r/M — 2). The field asymptotics 
at the horizon r — > 2M corresponds to an infalling wave 0(t + r#). The boundary conditions 
far away from the black hole follow from the assumption that there the cosmological field 
is spatially uniform and is free to roll down the field potential. We take 4>{t,r) = (pooif) at 
r ^> 2M, where the time dependent function 0oo(^) is the background evolution of the field 
in the absence of the black hole. The function 4>oo(t) depends on the scalar field potential 

v{<t>). 

Equation (jT^) is well studied for the case of a massive scalar field V(<f)) = ^m 2 (f) 2 , where 
it can be reduced to a linear Schrodinger-type ordinary differential equation with an effective 
potential U e s(r). 

However, in general the equation ( |T0"D is non-linear and the behavior of its solution is 
very different from the case of a massive scalar field; compare Figures | and |. We treat 
equation (|10"D as a partial differential equation and solve it numerically for several examples 
of runaway potentials. We also consider the case of a very light massive field and compare 
it with the case of a heavy massive field. 

The sequence of radial profiles <j>(t,r) for several moments t is shown in Figure |] for a 
runaway potential and in Figure ^ for a very light massive field, respectively. 

We found that r) in these cases with a very high accuracy can be approximated by 
a simple formula 



(f)(t,r) = 4> 



2M\ 

t + 2M\n (1 

r / 



(11) 



This result is rather insensitive to the form of the runaway potential. We will call the solution 
(|TT|) the delayed field approximation. This is because the field profile is merely delayed, not 
frozen, on the horizon. In the Eddington-Finkelstein coordinates (v , r) the incoming light 
geodesic is v = t + r* . The solution fllTD near the horizon is ~ 4>{v — vq) so the field crosses 
the horizon at the speed of light, but is delayed by a small amount vq = (2 In 2 — 1)M 
compared to the null characteristic. Far away from the black hole, at large r 3> 2M, a 
long-range tail is formed, which can be described by the formula 



<p(t,r) w 0ooOO 



(2Mf 



(12) 



For a light massive field when the equation (|10"D is linear, the tail ([12]) can be related to 
scattering on the effective potential U e s(r) with a sharp peak at ~ 3M. 

Our approximation (0), verified by the numerical solutions, generalizes the similar an- 
alytic result of for the case of a free scalar field (V = 0), where (fioo(t) = t. The 
accuracy of the delayed field approximation can be estimated by substituting the expression 
(p|) into the differential equation (|T0|) and calculating the residual term which is of order 
O f^r-(l — ))■ The residual is bounded and localized at around 3M; this is why the 
approximation (O) works. 

Equipped with equation (|TTD, we can calculate the energy flux of the scalar field through 
the horizon and the accretion rate onto the black hole M = ^4f. We can calculate M 

at 
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FIG. 2. Field profiles <p(t, r) outside a black hole for the runaway potential V(cfi) = l/(/> 2 at two 
subsequent time moments; horizontal axis is a tortoise coordinate r*. Red (solid) curves are the 
numerical solutions, and green (pale) curves are obtained with the delayed field approximation. 
The horizontal line corresponds to the homogeneous initial conditions. 
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FIG. 3. Field profiles for the potential V((f>) = ~m (f) with a very small mass m at two 
subsequent time moments. Red (solid) curves are the numerical solutions, and green (pale) curves 
are the delayed field approximations. 




FIG. 4. A snapshot of the field profile for the potential V{4>) = km <jr with a large mass m 
when the field oscillates. Red (solid) curve is the numerical solution, and green (pale) curve is the 
delayed field approximation. g 



in different coordinates. In Schwarzschild coordinates, the spatial gradient of develops 
just outside the horizon, and the energy flux is defined by T r * = 4>j<t>, r ~ 0^. In these 
coordinates, an observer at asymptotic infinity does not actually see the energy flux go 
through the horizon. What she/he sees is the energy from accretion accumulating just 
outside the Schwarzschild radius, increasing the black hole mass (defined in the ADM sense). 
In the Eddington-Finkelstein coordinates, the energy flux through the horizon is given by 
T vv = <i>l. 

Independent of the choice of coordinates, the rate of accretion of the scalar field by a 
black hole and the rate of increase of its mass in the regime of the delayed field approximation 
is equal to 

M = 4tt(2GM) 2 0^ (13) 

It is easy to check that this formula is equivalent to the first law of black hole thermodynamics 
8E = TdS. We can rewrite (|13|) in the form (jg^J = —AirGcj) 2 ^, which is readily similar 
to the corresponding formula (|6]) for the slow roll scalar field in quasi-de Sitter geometry. 

If we turn to the case of a massive scalar field with V(0) = |m 2 2 , when the background 
scalar oscillates around the minimum, we find that the delayed field approximation fllTD 
quickly fails, see Figure |], and formula fll3p is not valid. In this case the scalar field can 
be described in the WKB approximation. The accretion of the massive scalar field looks 
like accretion of free massive quasi-particles of mass m. Notice that this is different from 
a massive oscillating scalar field in the de Sitter background, where the flux of the heavy 
scalar field through the horizon remains the same as for the light slow roll field, ~ <f) 2 . 

In Appendix A, we briefly discuss some astrophysical applications of the results of this 
section. 



IV. PERTURBED METRIC AND HORIZON THERMODYNAMICS 

Classical perturbations in the de Sitter geometry die out in accordance with the "no- 
hair" theorem. Therefore, as in the black hole case, there are well-defined concepts of the 
entropy S and temperature T associated with the pure de Sitter horizon [13| . 



The spherical symmetry of the spacetime naturally selects a foliation by surfaces of the 
constant physical radius p = const, with normal vector V M p. If the spacetime is static (i.e 
admits a timelike Killing vector field fc M with vanishing differential invariant k li;u kpe IJ ' l '" a \ do 
not mix with static form of the metric), as de Sitter spacetime is for example, the vector 
£ M = —e^VvP orthogonal to V M p is also a Killing vector and a null generator of the event 
horizon (which coincides with the apparent and Killing horizons in the static case). Therefore 
one can define the surface gravity at the horizon k by ^"^u = and identify k = H and 
the temperature in the usual way T = The area entropy is equal to 1/4 of the apparent 
horizon area in the Planck units, S = 

For the quasi-de Sitter geometry with the homogeneous scalar field <p(t) slowly rolling 
towards its minimum, when the horizon radius 1/H(t) is adiabatically changing with time, 
we still can introduce an entropy S = ^ j^h^ and temperature T = associated with the 
apparent horizon. 
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In inflationary cosmology, classical perturbations which could exist prior to inflation die 
out. However, one of the striking features of the quasi-de Sitter inflationary stage is that the 
quantum fluctuations of the light (minimally coupled to gravity) scalar field 8(f) are unstable 
and inevitably generate scalar metric perturbations $(t, x) so that the quasi-de Sitter 

metric ([I]) becomes perturbed 

ds 2 = -(l + 2$)dt 2 + (l-2$)e 2 J dtH (dr 2 + r 2 dtf) . (14) 

If the spacetime is not static (timelike Killing vector does not exist), as in the inflationary 
geometry (|14|), the situation with the horizon thermodynamics is not quite so simple. The 
event horizon and apparent horizon are in general different surfaces, and the Killing horizon 
does not exist, so the notion of surface gravity is ill-defined. 

Fortunately, in the case of spherical symmetry one can define the mass (energy) M inside 
the spherical region of the physical radius p by 1 — = (Vp) 2 . The mass function M 
obeys the mass formula |TjJ which follows from the Einstein equations 

M, a = Arrp 2 (T ab -Tg ab )p'\ (15) 

where a, b stands for (t, r) only. We derive the mass formula ( [To] ) in Appendix B and give 
there other definitions to be used in this Section. 

We will argue now that the mass formula flip] ) applied at the apparent horizon can 
be interpreted as a first law of thermodynamics 5E = TdS for an arbitrary spherically 
symmetric metric. In the next section, we will apply this result to the perturbed inflationary 
metric (|14|), which as we will show can be well treated as (locally) spherically symmetric. 

Let us look at the apparent horizon, as it is locally defined and much easier to find than 
an event horizon. In spherical symmetry, the position of the apparent horizon is given by 
/ = (Vp) 2 = 0. The vector V a / is normal to the apparent horizon, which is the surface 
/ = 0, while the orthogonal vector ( a = £ ab Vb/ is tangent to it. Unlike in the static case, 
these vectors are not necessary null. The change of the mass function M along the apparent 
horizon is trivially related to the change of the radius of apparent horizon, 

G( a M, a = K P Cp, a = ± C(vrp 2 );a , (16) 

which allows its interpretation as the first law of thermodynamics 

rp dS 

^= < i7) 

The parameter k is defined in the Appendix B and generalizes the surface gravity. To relate 
the heat flow term with the stress-energy tensor, we use the mass formula (15) which gives 



CM, a = ^p 2 {T b a -Tg b a )C P , b . (18) 

The above expression can be rewritten in a more convenient form if one realizes that ( a M. a = 
-2Ki a M, a 
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- c (scL= 2 - 4r '^ (i9) 

This formula relates the change of the apparent horizon area to the (spherically symmetric) 
flux of matter through the p = const surface passing through the point on the horizon at 
which M is evaluated. 



V. FLUCTUATIONS FROM INFLATION AND WIGGLES OF THE HORIZON 

AREA 

In the time-dependent metric ([[]) of the quasi-de Sitter geometry, quantum fluctuations 
of the inflaton scalar field can be expanded with respect to the eigenmodes 

6<f>(t, X) = J dk {®klm 0kJm + atlm $Mm) > ( 20 ) 
Im 

where a^im and djy m are annihilation and creation operators; 4>ki m = 4>k(t) 2kji(kr) Yi m (Q) are 
the eigenmodes. The r— dependent eigenmode factor is given in terms of the spherical Bessel 
functions ji(kr), Yi m (Q) are the spherical harmonics, and the time-dependent eigenmode 
factor is given in terms of Hankel functions (j>k{t) = ^ Hrf/ 2 TC^(kr]), where rj is the 
conformal time r\ = —jjZ~ Ht - The Bunch-Davies vacuum corresponds to the absence of 
particles in the past dfc/ m |0) = 0. This corresponds to the positive-frequency asymptotic 
fc (t) ~ -^=e~ lhv in the far past rj — > — oo. Quantum fluctuations of the inflaton field 
are unstable and turn into long-wavelength classical inhomogeneities. Indeed, the modes 
which have physical wavelengths smaller than the horizon 1/H, have krj ^> 1 at the horizon 
and quickly oscillate without producing physical effects. In contrast, the modes with the 
wavelengths which exceed the horizon size kt] = k < 1 cease to oscillate, freeze out, and 
look like a classical scalar field with amplitude 0fc(t) ~ ^3/2 ■ 

We need to know the form of the eigenmodes at the apparent horizon where r = 77. 
Again, the time-oscillating modes with small physical wavelengths kr = kr] ^> 1 are spatially 
oscillating functions ji(kr). Large wavelengths modes with kr = kr] 1 have asymptotics 
ji(kr) ~ (kr) 1 . Notably, only the s-wave with / = survives in this asymptotic. Therefore, in 
our discussion of the wiggles of the horizon area due to metric fluctuations, we can consider 
only spherically symmetric perturbations. 

Fluctuations 5<j) generate scalar metric perturbations which are represented in flliD by 
scalar $(t,r). Now we are going to relate 5(f) and $ using the apparent horizon thermody- 
namics. 

In the perturbed spacetime (|TJp, the physical radius is p = (1 — Q)ef dtH r. Apparent 
horizon is defined by (Vp) 2 = 0. From here we find for the apparent horizon 

* = ir-'(i + . + |-£). (21) 

The entropy of the apparent horizon and the temperature are identified with geometrical 
quantities as 
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where k is denned by equation (|3"4]) of Appendix B. Evaluating the equation ( fL9f) relating 



the change in the horizon area to the energy flux at the position of the apparent horizon in 
the perturbed spacetime Ql4]) , and keeping only the linear terms with respect to $, 5(f), and 
their derivatives, one finds the following expression for the perturbations 

$ - - $ r + \ $ rr + H§ = 8nG <p(S(p-- Sip,) , (23) 



a a \ a 



valid at the position of the apparent horizon. 

Even though we know already that the first law of thermodynamics fll7D originated from 
the Einstein equations in the form of (|T5|), and hence the above equation would follow from 
the linearized Einstein equations, let us explicitly demonstrate this fact. The usual equations 
for metric perturbations from scalar field fluctuations are |2E 

I ' ■ V| ^=(^V , W=4*Gto>, (24) 



y 47rGa 2 ip 2 / 

where the for spherically symmetric $, V|$ = <9 rr $ + ^<9 r $. Combining both equations by 
substituting them into an identity 2ip5(p = (ft 2 + (<p5<p), we find that 

$--$ r + \§ rr + H$ + ^-(--a)$ r = 8nG <h ( 5(2> 5ip r ) . (25) 

a ' a 2 ' a 2 \r ) \ a / 

This equation holds everywhere in the perturbed spacetime. When evaluated at the apparent 
horizon, the last term in the left-hand side vanishes and we are left exactly with the equation 

(ED- 

Consider the area of the apparent horizon A = 47rp^. For frozen long- wavelength fluctu- 



ations, we can drop derivatives in formula (|2T| ) for p^. Then the leading contribution of the 
frozen classical fluctuations to area is 

A=— (1 + 2$) . (26) 

Thus, the quantum fluctuations of the inflaton field generate wiggles of the quasi-de Sitter 
horizon ctrGclj clS described by fl2"6| ) and sketched in Figure [5[ 

Locally (within the Hubble patch), almost spherical metric perturbation $ can be viewed 
as the variation of the Hubble parameter H PJl2]|. Indeed, transform the metric ( |14|) to 



synchronous coordinates i = t + J* 1 dtQ, x, = Xi + J — J Vj$. Then we have ds 2 = —dt 2 + 
g 2 J dtii ^ w ]j ere ^ ne new Hubble parameter is H — (1 — $)iJ. As expected, the area of the 
horizon is the same 

4:71 4:71 

A = jp( 1 + 2 * ,) = T 2 - (27) 

Therefore an observer within a Hubble patch sees the local Hubble parameter H simultane- 
ously slow rolling and wiggling due to the instability of quantum fluctuations of the inflaton 
field, in the spirit of stochastic picture of inflation. 
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VP 





FIG. 5. A sketch of the apparent horizon of radius p. Left: the horizon in quasi-de Sitter 
geometry where p = 4j very slowly increases with time. Right: the same but with small metric 
perturbations where p = -^(1 + $) is wiggling with time due to the randomly varying <£. 



We have to note that there is another, much smaller correction to the geometrical entropy 
due to the back reaction of quantum field effects in curved background. Indeed, quantum 
gravity effects give one-loop corrections of order of R 2 ~ H 4 in the right hand side of the 



Einstein equations, which slightly shift the Hubble parameter [22] and the horizon area by a 



factor l + c-S, where c is some numerical coefficient. However, this leads to a tiny correction 



in the area entropy. 

Let us discuss the formula 



C 



y3/2 



1 



-p v ,4> 

of the inflaton field 



where C 



). Metric fluctuations $(£, r) ~ 5 
(37r) 3 / 2 is a numerical coefficient. 



• have a flat spectrum 
Classical fluctuations 



can be treated as a random (Wiener) process with zero vacuum 



expectation value, but with the increasing dispersion (5 (ft 2 ) ~ f aH * = f^t. Suppose 
for a moment that the background value of H(t) is changing very slowly. Then, scalar 
metric fluctuations r) have zero mean value but their dispersion is increasing with 
time linearly, ($ 2 ) = C 2 J aH f ~ t. Thus, the area of the apparent horizon looks 

P ,4> 

like a random (Wiener) process with the mean value A = jMvs and with the dispersion 
($ 2 ) linearly increasing with time. Imagine an ensemble of inflationary Hubble patches with 
apparent horizon areas given by the formula (|26|). Then their area is a statistical value which 

(A-A) 2 

obeys the Gaussian distribution P(A, t) ~ e 2 <* 2 > . Now, if we take into account variation 
of background values like H(t) and V(4>(t)) with time, time-dependence and statistical 
properties of the horiazon area given by ( p6|) or equivalently (p7|) will be more complicated. 
We will discuss this issue below. 



VI. DISCUSSION: ENTROPY AND COSMOLOGICAL FLUCTUATIONS 

In this paper we calculated the apparent horizon area and geometrical entropy of quasi- 
de Sitter geometry, which describes an inflationary stage, and the energy flux of the scalar 
fields through the apparent horizon. Issues related to de Sitter thermodynamics are com- 
monly considered in the static de Sitter coordinates, where an observer at the origin is 
surrounded by the event horizon and detects a thermal flux of temperature T — 
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We work with the geometrical entropy S of the apparent horizon in the time- dependent 
planar de- Sitter coordinates adopted in inflationary cosmology, which admit a simple gener- 
alization to the quasi-de Sitter geometry and, most importantly, admit a clear interpretation 
of fluctuations generated from inflation. 

The slow roll of the inflaton field leads to slow change in the horizon radius. Assuming 
that the energy flux of the rolling scalar field through the horizon changes the geometrical 
entropy 8E = TdS, we reproduce the Einstein equation (0) which relates H and 2 . Other 
background scalars \ which are subdominant during inflation give similar contributions to 
the energy flux y 2 , independent of their potentials. 

This change of the quasi-de Sitter horizon radius due to the rolling scalar field is very 
similar to the increase of the mass of a black hole due to the accretion of a background scalar 
field rolling towards the minimum of its potential M ~ 2 . This type of accretion, which can 
be described analytically with the delayed field approximation, can be realized for runaway 
potentials of the background scalar field or for very light massive fields. We present this 
material here mainly to illustrate similar calculations for rolling scalars in quasi-de Sitter 
geometry, as the astrophysical effects of accretion of a rolling scalar onto a black hole that 
we considered (say quintessence and astrophysical-size black holes) are negligibly small. 
Oscillating heavy scalar field accretes onto black hole very differently as quasi-particles. 

Inflationary, quasi-de Sitter stage erases classical inhomogeneities which could exist prior 
to it. However, quantum fluctuations of the inflaton field (as well as other light scalars 
minimally coupled to gravity) are unstable and produce long-wavelength fluctuations of the 
inflaton field 50 which behave as classical inhomogeneities at scales larger than the Hubble 
patch of size if -1 . Fluctuations 5<p generate long- wavelength scalar metric perturbations $. 
Geometrical quantities, like an area of the apparent horizon, acquire corrections due to the 
scalar metric perturbations. We calculate the energy flux of the inhomogeneous scalar field 
<j)(t) + S(j)(t, x) through the apparent horizon and the change in the apparent horizon area 
of the perturbed metric multiplied by its (geometrical) temperature TdS. Again, equating 
5E = TdS, we show that this thermodynamics relation is compatible with the linearized 
Einstein equations which relate $ and 5(j). 

Thus, as long as the Einstein equations hold, generation of the inflaton fluctuations is in 
full agreement with the variation of the entropy of the quasi-de Sitter horizon. Therefore, 
in the picture of rolling inflaton with quantum fluctuations generated with time, we did 
not find that quantum fluctuations may violate holographic bound during inflation. (Notice 
that the fluctuations from inflation are described by squeezed states, which do not carry 
entropy p3|j24]] . In simple terms, locally the effect of fluctuations is just a wiggling of the 



local Hubble parameter). 

It was suggested in |1| that in the "hot tin can" picture the entropy of fluctuations may 
violate holographic bound during inflation.^ The issue of inflaton fluctuations in the "hot 
tin can" picture is not clear to us. In theory of inflation, scalar field fluctuations are usually 



x If frozen fluctuations from inflation were carrying large entropy, then a number N of free light 
scalars produced during inflation would have ./V times bigger entropy and UV cutoff proposed in 
[16] would depend on N. 
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considered in the time dependent metric (|l|), as it was described in Section [V[ Bunch-Davies 
vacuum corresponds to the thermal state in the static coordinates. On the other hand, in 
the static de Sitter coordinates, discussion of the quantum field theory is usually restricted 
to thermal radiation from the horizon, usually in terms of the detector response. Thermal 
radiation associated with the de Sitter horizon is related to any free quantum field, scalar 
fields with mass m and any coupling to gravity, vector fields etc., the difference will be 
only in the radial dependence of the wave functions ("gray-body factor"). On the other 
hand, instability of quantum fluctuations from inflation occurs only for very light minimally 
coupled scalar fields. It will be interesting to understand what is relation of the quantum 
field theory in the "hot tin can" picture, and the instability of fluctuations from inflation. 
For this it will be essential to work with regularized VEV < 5(p 2 >. 

As we already mentioned, in the picture of rolling scalar field with accumulating inflaton 
fluctuations, which we adopt in the paper, holographic bound and quantum fluctuations 
from inflation are compatible. 

A lesson from our calculations is that in the quasi-de Sitter geometry with slowly rolling 
scalar field, the horizon area, or geometrical entropy, are perturbed by the scalar metric 



fluctuations according to the formulas (p6[) , (|27|). This means that the area of the horizon 
for an ensemble of the different Hubble patches is not the same but is a statistical variable by 
itself. For small metric perturbations (of order of $ ~ its mean value is A = j^, where 
H{t) is slowly decreasing background value, but dispersion around the mean value is defined 
by ($ 2 ), which growth with time. Equivalently, we can talk about statistical properties of 
the local Hubble values H @Jl^]. Notably, this picture is converging with the stochastic 




description of inflation in terms of the probability distribution of the inflaton field P,|T!! 

Probability to have the value of inflaton field (f> in quasi-static regime is P{4>) ~ exp , srf> } 

and can be interpreted in terms of entropy S = logP(0) |p5| . Remarkably, this entropy is 
identical to the geometrical entropy During inflation and rolling of <fi the Wiener process 
of accumulating fluctuations 5<p (or similarly perturbations of $) changes distribution of 
P(0). Distribution of the horizon areas of different Hubble patches, defined by the local 
values of the Hubble parameters, depends on the background slow roll regime and the regime 
of accumulation of fluctuations, both of which depend on the model of inflation []9|, [T0|JT^] . It 
would be interesting to understand further the correspondence between stochastic approach 
to inflation and geometrical entropy (PB"D, (^7J). 

So far we discussed small metric fluctuations $ or small local variations of H. However, 
in the chaotic inflationary scenario for large enough values of </>, variations of H due to the 



quantum jumps of S(p can be large and lead to the self- reproducing inflationary universe [11 



We would like to draw attention to this regime (which is still below the Planck energy den- 
sity) and to note that "adiabatic" geometrical thermodynamics which we considered in this 
paper is not applicable here. Indeed, consider the Hubble patch where H is increasing due 
to the quantum jumps. Increase of H is not compatible with the classical Einstein equation 
Consequently, quantum jumps are not compatible with the horizon thermodynamics 
(of a single Hubble patch) since 5E = TdS does not hold either. Geometrical entropy of the 
local Hubble patch is decreasing. However, we have to take into account the entropy of all 
Hubble patches. Self-reproduction of inflating regions looks like a chain reaction, which is 



described by the branching diffusion process [26 
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APPENDIX A: ROLLING COSMIC SCALARS AND BLACK HOLES 

Here we briefly discuss astrophysical applications of the results of Section [TTI[ One of the 
most important examples of slowly rolling scalar field is a homogeneous cosmological scalar 
field, whose evolution is given by the field equation + 3H<f> + V^, = with the Hubble 
friction. This type of fields appears, for instance, in the cosmological models of quintessence. 
To dominate at the present cosmological stage while avoiding gravitational clustering, the 
value of (j> should be of order of the Planck mass M P , and its effective mass = very 
small, m^j < H, where H(t) in the present day Hubble parameter. 

Quintessence may interact gravitationally with black holes of different masses, ranging 
from tiny primordial black holes to astrophysical black holes of solar masses to supermassive 
black holes. Accretion of quintessence onto a black hole is given by the formula (|T3|). Scalar 
field velocity 0(t)oo is defined by its equation of motion; however, the kinetic energy <fi 2 
cannot exceed the total energy density of the universe MpH 2 . Substituting this upper limit 
in (|13D, we estimate the rate of the black hole growth due to the accretion of the rolling 
cosmological scalar field 

M ~ H 2 M 2 ~ (j^f , (28) 

where r g = 2M is the gravitational radius of the black hole, H~ x is the size of the universe, 
and we switched to Planck units Mp = 1. Thus, accretion of cosmic scalars including 
quintessence is absolutely negligible for all types of astrophysical black holes. For primordial 
black holes at the moment of formation, when the ratio \^=t) is of order of unity, the 
study of Section |TJ of rolling scalar field accretion is not applicable directly, as the black 
hole spacetime would be significantly different from Schwarzschild in this case. Even if 
one believes the equation flT3| ) in this regime, it would seem unlikely that the accretion of 
quintessence can cause explosive growth of primordial black holes (contrary to some claims 
in the literature). For this to occur, the seed primordial black hole would have to be much 
bigger than the cosmological horizon if the quintessence is subdominant. 

It is interesting to take into account the back reaction of an accretion on the evolution of 
the cosmic scalar field itself. From the energy balance we obtain a correction to the Hubble 
friction term due to the accretion of rolling scalar onto black holes 

+ (3H + w 9 2 )0 + V'(<t>) = , (29) 

where ubk is the spatial density of black holes. This correction is proportional to the filling 
factor of black holes, i.e. the fraction of volume occupied by black holes. The last is tiny so 
that the effect of the interaction of a rolling scalar field with black holes is negligible. 
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Rolling scalar field interacting with a black hole also emerges in the context of the Brans- 
Dicke or more general scalar-tensor theories |20|,19 . 



APPENDIX B: MASS FORMULA IN SPHERICAL SYMMETRY 

A general spherical spacetime is described by the metric 

ds 2 = g ab dx a dx b + p 2 dtt 2 , (30) 

where g ab is the metric on a two-manifold with coordinates (t, r), and p is the physical 
radius of spherical slices. The field dynamics is given by the Einstein-Hilbert action, which 
in spherical symmetry can be dimensionally reduced to yield 

S = J V=7 d 2 x [p 2 i?[ 7 ] + 2(Vp) 2 + 2] + p 2 £ mat t cr } • (31) 

The Einstein equations for spherically symmetric spacetime fl30l) follow from the above action 
by varying it with respect to the two-metric g a b and the two-scalar p. In particular, the {ah} 
component of the Einstein equations is 

2p- 1 [g ab Up - p. ab ] + p~ 2 g ab [{Vp) 2 - 1] = 87rGT afe . (32) 

By subtracting the (two-dimensional) trace, this equation can be written in a more conve- 
nient form 

P-ab - ig ab = -AixG p(T ab - Tg ab ), (33) 

where we have defined 

1 

2p 



K = -[l-(Vp) 2 ]. (34) 



In the case of Schwarzschild black hole, k coincides with the surface gravity at the horizon. 
We will continue to use the same notation since k defined in (|34|) will enter the thermody- 
namics relation in the more general case. 

In spherical symmetry, it is possible to define a local mass function M(x a ) by 

1 - ^ = (VP) 2 ^ / (35, 

The change in mass is related to the flux of matter given by the stress-energy tensor T ab . 
To see this, let us take the derivative of M — | (1 — /) 

GM. a = -|(/ ;a - 2Kp. a ) = -p(p- ab - K9ab)p' b (36) 

and use the Einstein equation (|33|) we derived above to obtain 

M, a = 47rp 2 (T ab -Tg ab )p' b . (37) 



This is the differential mass formula in a spherically symmetric spacetime WE . 
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